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Abstract 

The metallicity distribution in the intracluster medium of the NGC 5044 group was studied up to 
0.3ri8o using the XIS instrument on board the Suzaku satellite. Abundances of O, Mg, Si, S, and Fe were 
measured with high accuracy. The region within a radius of 0.05r 18 o from the center shows approximately 
solar abundances of Mg, Si, S, and Fe, while the O/Fe ratio is about 0.5-0.6 in solar units. In the outer 
region, the Fe abundance gradually drops to 0.3 solar. Radial abundance profiles of Mg, Si and S are similar 
to that of Fe, while that of O seems to be flatter. At r > 0.05rigo, the mass density profile of O differs from 
that of Fe, showing a shoulder-like structure that traces the luminosity density profile of galaxies. The 
mass-to-light ratios for O and Fe in NGC 5044 are one of the largest among groups of galaxies, but they are 
still smaller than those in rich clusters. These abundance features probably reflect the metal enrichment 
history of this relaxed group hosting a giant elliptical galaxy in the center. 

Key words: galaxies:abundances — clusters of galaxies: intracluster medium — groups individual 
(NGC 5044 group) 
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1. Introduction 

Groups and clusters of galaxies include most galaxies in 
the nearby universe. (Mulchaey 2000). These systems also 
represent building blocks of rich clusters and are the best 
laboratories for the study of thermal and chemical history 
governed by baryons. An important clue in studying the 
evolution of galaxies is the abundances of different ele- 
ments in the hot X-ray emitting gas in groups or clusters 
of galaxies: namely the intracluster medium (ICM). The 
metals in the ICM have been synthesized by supernova 
(SN) in galaxies. As a result, ratios of ICM metal mass 
to the total light from galaxies in clusters or groups, i.e., 
metal-mass-to-light ratios, are key parameters in investi- 
gating the chemical evolution of the ICM. 

Since Fe lines are prominent in X-ray spectra, the con- 
tent and distribution of Fe in the ICM have been studied 
in detail. The ASCA satellite first enabled us to measure 
the distribution of Fe in the ICM (e.g., Fukazawa et al. 
2000; Finoguenov et al. 1999). Poor clusters and groups 
of galaxies exhibit different properties from richer systems 
in that the derived iron-mass-to light ratios (IMLR) are 
systematically smaller (Makishima et al. 2001). 

Rasmussen & Ponman (2007) reported the temperature 
and abundance profiles of 15 nearby groups of galaxies ob- 
served by Chandra. They showed that Fe abundances of 
the groups declined from about 1 solar at the center to 
about 0.1 solar at r 50 o, which is about 0.54 ngo- Si shows 
less systematic radial variation, indicating SN II contribu- 
tion increases with radius. Finoguenov et al. (2007) also 



reported the metal distributions of groups of galaxies ob- 
served by XMM-Newton. They also found that metallicity 
of some groups drop to 0.1 solar at ~ r 500 . 

O and Mg are synthesized predominantly in SN II, 
while Fe and Si are synthesized in both SN la and SN II. 
Abundance measurements spanning the range of species 
from O to Fe are therefore required for unambiguous de- 
termination of the formation history of massive stars. 
XMM-Newton provided the means to constrain O and 
Mg abundances in some systems (e.g., Matsushita et al. 
2003,Tamura et al. 2003, Matsushita et al. 2007b), how- 
ever, reliable results have been obtained only for the cen- 
tral regions of very bright clusters or groups of galaxies 
dominated by cD galaxies. 

Suzaku (Mitsuda et al. 2007) is the fifth Japanese 
X-ray astronomy satellite. The X-ray Imaging 
Spectrometer(XIS) instrument (Koyamaet al. 2007) offers 
an improved line spread function due to a very small low- 
pulse-height tail in the energy range below 1 keV coupled 
with a very low background. Therefore, especially for the 
regions of low surface brightness or equivalent width, the 
XIS provides better sensitivity for O lines. The instrumen- 
tal Al line of the MOS detectors on XMM-Newton causes 
problems in measuring the Mg abundance in somewhat 
fainter systems. 

With the Suzaku satellite, the oxygen-mass-to-light ra- 
tio (OMLR) as well as the IMLR of several clusters of 
galaxies and several groups of galaxies were measured up 
to 0.2-0.3ri 80 (Matsushita et al. 2007a, Tokoi et al. 2008, 
Sato et al. 2007a, Sato et al. 2008a, Sato et al. 2008b, 
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Sato et al. 2008c) . The OMLR and IMLR increase with 
the radius, and become relatively flat beyond 0.1 riso- 
The Fornax cluster, which is the nearest poor cluster with 
asymmetric X-ray morphology, shows the smallest IMLR 
and OMLR up to 0.13 r 180 (Matsushita et al. 2007a). 

The NGC 5044 group is a nearby group of galaxies 
whose redshift is z = 0.009020 from the NASA/IPAC 
Extragalactic Database (NED). The group has a giant el- 
liptical galaxy NGC 5044 at the center. The X-ray emis- 
sion of the group shows nearly symmetric spatial distribu- 
tion (David et al. 1994). The ICM properties in this group 
have been studied well using XMM-Newton and Chandra 
by Buote et al. (2003a) and Buote et al. (2003b), and the 
metal distribution in the outer regions was studied using 
XMM-Newton (Buote et al. 2004). They revealed that 
the Fe abundance of the ICM drops from about 1 solar 
within r « 50 kpc to ~ 0.4 solar near r ~ 100 kpc using 
the solar abundances by Grevesse & Sauval (1998) (Buote 
et al. 2003b). In the outer region, at ~ 176-352 kpc, 
the iron abundance is determined to be about 0.15 so- 
lar using the same solar abundances (Buote et al. 2004). 
Tamura et al. (2003) showed O abundance to be 0.25 solar 
within the central 10 — 20 kpc using the Reflection Grating 
Spectrometers (RGS) on XMM-Newton. The mass distri- 
bution has been studied using Chandra and ROSAT data 
(David et al. 1994, Betoya-Nonesa et al. 2006). 

We use the Hubble constant H = 70 km/s/Mpc. The 
distance to the NGC 5044 group is D L = 38.9 Mpc, and 
1' corresponds to 11.1 kpc. The virial radius, riso = 
1.95 (ff /100)- 1 v /fe{T)7r0 keV Mpc (Markevitch et al. 
1998, Evrard et al. 1996), is about 880 kpc for the av- 
erage temperature k(T) = 1.0 keV. We use the new 
abundance table from Lodders et al. (2003) in this pa- 
per. Abundances of O and Fe are about 1.7 times and 
1.6 times higher than those of Anders & Grevesse (1989), 
respectively. Unless otherwise specified, errors are quoted 
at 90% confidence for the single parameter of interest. 

2. Observations 

Suzaku performed three pointing observations of the 
NGC 5044 group in July 2006. The observational log 
is shown in Table 1. The first observation (hereafter, 
the central field) had a pointing direction toward the cD 
galaxy, NGC 5044. The second and third observations 
were centered 15' north and 15' east of NGC 5044, respec- 
tively (hereafter, the north and east fields, respectively). 
Figure 1 shows a 0.5-4.0 keV image for the three fields. 
The observed region covers a distance of about 24 arcmin, 
or - 267 kpc, from NGC 5044. 

We also analyzed the MOS and PN data from the 
XMM-Newton observation of NGC 5044. Details of the 
data reduction are given in Nagino et al. (2008). 

3. Spectral Analysis 

The XIS was operated in the nominal mode. We in- 
cluded the data formats of both 5x5 and 3x3 editing 
modes in our analysis using XSELECT(Ver . 2.3). The 
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Fig. 1. XIS image (0.5-4.0 keV) of the NGC 5044 Group. 
The NXB (non X-ray background) and CXB (cosmic X-ray 
background) were subtracted, and the difference in expo- 
sure times was corrected. The image was smoothed with a 
Gaussian of a = 16 pixels. 

analysis was performed using HEAsof t (Ver . 6.0.6) and 
XSPEC (Ver . 1 1 . 3 . 2o) . After applying the standard data 
selection criteria, the exposure times of center, north, and 
east fields were 19.7, 54.6, and 62.4 ks, respectively. 

The response of the XRT and XIS was calculated using 
the xissimarf gen ancillary response file (ARF) generator 
(Ishisaki et al. 2007) and xisrmf gen response matrix file 
(RMF) generator, version 2006-10-26. Slight degrada- 
tion of the energy resolution was considered in the RMF, 
and decrease in the low-energy transmission of the XIS op- 
tical blocking filter (OBF) was included in the ARF. The 
ARF response was calculated assuming a surface bright- 
ness profile, S(r), based on the analysis of XMM-Newton 
data by Nagino et al. (2008). 

For the central field, we analyzed the spectra extracted 
from four annular regions of 0-2', 2-4', 4-6', and 6-9' cen- 
tered on NGC 5044. For the north and east fields, we used 
the spectra extracted from two annular regions of r < 15' 
and r > 15' centered on the galaxy. Each spectrum was 
binned to observe details in metal lines, especially O lines. 
Each spectral bin contained 50 or more counts. 

The non X-ray background (NXB) was subtracted from 
the spectra using the database of night Earth observations 
with Suzaku for the same detector area and with the same 
distribution of COR (Tawa et al. 2007). We produced 
the spectra for an energy range 0.4-7.0 keV, since above 
7.0 keV, subtraction of background lines is difficult. The 
spectra of XIS detectors (XIS0, XIS1, XIS2, XIS3) were 
fitted simultaneously with a model consisting of the ICM 
emission with additional power-law and the Galactic com- 
ponents. 



No. ] Suzaku observations of the NGC 5044 group 3 

Table 1. Suzaku observations of the NGC 5044 group 

Fields Seq. No. (RA, Dec) in J2000 Date of obs. Exp, time 

Center 801046010 (13 h 15 m 24. s 0, -16°23'07."8) 2006/07/02 19.7 ks 

North 801047010 (13 h 15 m 24?0, -16°08'07."8) 2006/07/03 54.6 ks 

East 801048010 (13 h 16 m 26?5, -16°23'07."8) 2006/07/04 62.4 ks 



3.1. The ICM component 

First, we assumed the ICM consisting of a single tem- 
perature vapec (Smith et al. 2001) model (hereafter, the 
IT model). The metal abundances of He, C, N, and Al 
were fixed to the solar values. We divided the other met- 
als into seven groups, O; Ne; Mg; Si; S, Ar and Ca; Fe; 
and Ni, and allowed them to vary. The spectra were sub- 
ject to a common interstellar absorption, Nn, fixed at the 
Galactic values, which are 5.03 x 10 20 , 4.88 x 10 20 , 5.21 x 
10 20 (atoms/cm 2 ) for the center, north, and east fields, 
respectively. 

We also applied a two-temperature (hereafter, 2T) 
vapec model for the ICM, where the metal abundances 
of the two components were assumed to have the same 
value. 

3.2. CXB and discrete sources 

The cosmic X-ray background (CXB) was modeled by 
a power-law spectrum with a photon index V = 1.4. The 
normalization was allowed to vary, which can also account 
for contributions from discrete sources in galaxies and a 
possible hard emission. 

3.3. The Galactic component 

The Galactic emission mainly arises from the local hot 
bubble (LHB) and the Milky- Way halo (MWH). We used 
two-temperature apec thermal spectra for this emission, 
and by combining it with the IT model for the ICM we 
fitted the observed spectra for the outermost (r > 15') 
regions in the north and east fields. The temperature 
and normalization of the two components in the Galactic 
emission were left free with metal abundance fixed to the 
solar level. The results of the spectral fits are shown in 
Table 2. The temperature values of the two apec models 
are consistent with the typical values for the LHB and 
MWH derived by XMM (Lumb et al. 2002). 

Since the errors for these temperatures are small, we 
fixed the temperatures of the LHB and MWH to 0.12 keV 
and 0.30 keV, respectively. Assuming that the Galactic 
component has the same surface brightness within the en- 
tire observed region, we calculated the normalization of 
the Galactic component in each region based on the sur- 
face brightness measured in the outermost region in the 
north field. The spectra of the other regions were con- 
sistently reproduced with the sum of the ICM model and 
this Galactic component. This feature supports the view 
that these apec components in the outermost region are 
dominated by the emission in our Galaxy. 



Table 2. Temperature of the Galactic components with the 
IT model for the ICM. 



Region kTum* HW Ratio of xVd.o.f. 

(arcmin/riso) (keV) (keV) Norm'' 

> 15/0.19 (N) O.Utooi - 30 -oo3 0.481^1 1113/895 

> 15/0.19 (E) 0.13±g;gi 0.29±g;°^ 0.471^3 996/895 



t Ratio of normalization iVormMWH/-'V orm LHB 

* LHB = Local Hot Bubble, MWH = Milky Way Halo 

4. Results 

4-1. Fitting results 

Table 3 summarizes the best-fit parameters of tempera- 
ture, x 2 values, degrees of freedom, and abundances. For 
the 2T fitting, the ratios of the vapec normalizations and 
flux between cool and hot ICM components are also sum- 
marized in Table 3. 

Within r = 0.050 ri 80 , or 4', the 2T model significantly 
improved the fit. For 0.050 r 180 < r < 0.19 r 180 (4' < r < 
15'), the 2T model provided a slightly better fit than the 
IT model, while at r > 0.19 riso (15') the \ 2 values derived 
for the IT and 2T models were nearly equal. 

In Figure 2, representative spectra obtained using the 
XIS-0 and XIS-1 instruments are shown. The spectra are 
well fitted with the 2T model, and the IT model provides 
reasonable fits in the outer region. 

4-2. Temperature profile of the ICM 

We derived the radial temperature profile based on the 
spectral fit for the annular regions from the group center 
to r ~ 0.3 ri 80 (Figure 3a). The fit with the IT model, 
within r < 0.1 riso, showed a temperature increase from 
0.8 keV to 1.1 keV. At r > 0.2 r 1S0 , it dropped slightly to 
1 keV. 

For the 2T model, the temperature of one component 
showed similar features to that for the IT fit. Within 
r < 0.1 ~ 0.2 r 180 , we required another component with 
kT ~ 1.5 keV. The normalization ratio for the cool and 
hot vapec components is shown in Figure 3 (b). The cool 
component is the strongest in the innermost region. In 
the outer region, the two components show similar nor- 
malizations. 

4-3. Abundance profiles of the ICM 

Radial profiles up to r < 0.3 r^o, or 260 kpc, of abun- 
dances of O, Ne, Mg, Si, S, Ar and Ca and Fe are sum- 
marized in Table 3 and Figure 4. 

On the whole, abundances are high at the center, and 
decrease toward the outer region. At the center, the IT 
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Fig. 2. Representative spectra from the XIS-0 (red) and the XIS-1 (black) instruments fitted with the IT (blue) or 2T (blue and 
light blue) models for the ICM including the power-law (orange) and the Galactic components (magenta and green). The lower 
panels show the fit residuals in units of cr. 
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Fig. 3. (a) Temperature profile of the ICM derived with the IT (gray) and the 2T (black) models for the center and the north 
fields (crosses) and for the cast field (diamonds), respectively. Dotted lines indicate shifts of the best-fit values when the NXB level 
is changed by ±10%. The dashed lines denote shifts when the thickness of the OBF contaminant is changed by ±10%. (b) The 
normalization ratio of the cool and hot components derived from the 2T model fit. 




Fig. 4. Abundance profiles of (a) O, (b) Ne, (c) Mg, (d) Si, (e) S and (f) Fe of the ICM derived from the IT (gray) and 2T (black) 
models for the center and north fields (crosses) and for east field (diamonds), respectively. The dotted lines indicate shifts of the 
best-fit values when the NXB level is changed by ±10%. The dashed lines denote shifts when the thickness of the OBF contaminant 
is changed by ±10%. 
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Table 3. The ICM temperature, 
the NGC 5044 group. Errors arc in the 90% confidence range of statistical errors, and do not include systematic errors. 



Fitting 
Region 
(arcmin/ri 8 o) 



IT model for the ICM 



2T model for the ICM 



kT 
(keV) 



XVd.oi. 



u. iyo_o .oo6 
0.961™! 

-+0.01 
-0.01 

1 06+ 02 

l-UO-0.02 



(keV) 



n 77q +().0UY 
u -< < y -o.oo7 

n 897+0-007 
U.OZ I _ .008 

96+ 05 

u - yo -0.12 

1 m+ 003 

-1.U1_q.08 



kT hot 
(keV) 



+0.10 
0.06 
+0.06 
-0.08 
h0.18 
-0.18 
1 74+0-28 
1 - ' 4 -0.21 



Ratio of 
Norm) 

Z - Oy -0.33 

°- 88 -0:09 
n 79 +0.30 
u - ' z -0.26 
n 77 +0.26 
u - ' ' -0.27 



Ratio of 
Flux* 



X 2 /d.o.f. 



0-2/0-0.025 
2-4/0.025-0.050 
4-6/0.050-0.076 
6-9/0.076-0.11 



1.05] 



i ' uo -0.02 

1 06+ 02 

1 - UD -0.02 

qo+ 04 

u - yu -0.05 
n Q7+0-05 



1574/899 
1441/899 
1194/899 
1027/899 



1.57 
1.41 
1.48 



3.88 
1.25 
0.90 
0.94 



+0.03 
0.04 
0.02 
0.05 
0.02 
-0.02 
0.09 
0.07 



+0.64 
-0.41 
-0.39 
-0.24 
-0.04 
-0.04 
+78.13 
1.57 



i-44±8:S 

0-89t° o f 2 

o nc+3.40 
J - UO (-3.05) 



1276/897 
1206/897 
1123/897 
967/897 



< 15/0.19 (N) 

< 15/0.19 (E) 

> 15/0.19 (N) 

> 15/0.19 (E) 



1047/899 
963/899 
1091/899 
970/899 



1.01 
1.02 
0.24 
0.91 



2.14 
1.72 
0.95 
1.77 



1.32 
1.61 
0.27 
2.95 



IT model for the ICM 



999.5/897 
927/897 
1112/897 
968/897 



Region 
(arcmin/rigo) 



O 

(solar) 



29 _FTrnr 

u - zy -0.09 

4Q+ 013 

u -^ y -0.12 

o.i5_8:iS 

OjOlgjf 

-0.15 
,+0.16 
3 -0.13 

<0.12 
<0.32 



Ne 
(solar) 



0.86™ 

n R q+0.29 

0.63_ 27 
7fi+°- 44 

U. IO_ 43 

-+u:33 

-0.34 
, +0.38 
-0.38 

0.59+ 029 



Mg 
(solar) 



Si 
(solar) 



-0.07 
-0.06 
0.06 



S 

(solar) 



-0.15 
-0.14 
0.12 



Fe 
(solar) 
0.65 +,UM 
0.62 

0-48-0.05 

n qr+0.04 
u -°°-0.04 

n 97+U.U4 
u.z^_ 04 

+0.04 



Ni 
(solar) 



FD733 
J -0.31 
9 oc+0.32 
z - o,J -0.30 

1 64+ 36 
1 -°^-0.32 

1 co+0.37 

1 - O J -0.32 

2+0.28 

-0.25 



0-2/0-0.025 
2-4/0.025-0.050 
4-6/0.050-0.076 
6-9/0.076-0.11 



0.94 
0.86 
0.69 
0.69 



-0.08 
5+0.09 
-0.08 
1+0-12 
J -0.11 
^+0.13 
-0.12 

o-eo±8:ii 

0.38l^ 

7+0.19 
-0.16 



0.70 

0-59_ . 06 

^J-'il -0 .07 
n o 7 +0.07 
U - J ' -0.07 
,+0.06 
-0.05 
+0.07 
-0.06 

1 5 +0 ' 07 

u - lo -0.07 
+0.11 
-0.10 



1.22 

0.95_ 011 

71+0 13 
u. 1 1_ 12 

52+ ' 14 

u - oz -0.13 
,+0.11 
-0.11 
+0.13 
-0.12 

S5+ ' 20 

3+0.30 
-0.27 



-0.04 
0.05 
0.05 

+0.05 



3.10] 



< 15/0.19 (N) 

< 15/0.19 (E) 

> 15/0.19 (N) 

> 15/0.19 (E) 



0.32] 
0.13 H 



0.65] 
0.91 H 



0.45 



0.30 
+0.48 
-0.45 



0.37] 



0.24] 
0.39 



0.12 



0.43] 
0.56 



0.43] 



0-34_ . 04 
22 +o.o6 



0.28 



-0.05 
+0.09 
-0.08 



2T model for the ICM 



0.98] 



1 or+0.33 

l.DO_ 29 

n qq+0- 36 

u - yy -0.35 
h0.65 

-0.55 



0.83 H 



Region 
(arcmin/ri 8 o) 



O 

(solar) 



Ne 
(solar) 



Mg 
(solar) 



Si 
(solar) 



S 

(solar) 



Fe 
(solar) 



Ni 
(solar) 



0-2/0-0.025 
2-4/0.025-0.050 
4-6/0.050-0.076 
6-9/0.076-0.11 


0- 58 i(x:L5 
76+ 21 

u -'°-0.19 

0.28_8;g 
S3+ 29 

u - oo -0.25 


1 77+0.43 
i -< '-0.35 

9 9/1+O.56 
z - z -0.55 

i-62l8:l 

l-80±°;?o 


1 ri+U.18 
1 - ol -0.16 
1 ci +0.20 
i - ol -0.18 

o.98_8;8 
o.98l8 : 2 g 


i.ool8;^ 
0-9218:11 
o.60l8;lo 

47+0-1 2 
u -*' -0.10 


1 01+0.19 

i - oi -0.17 

1 10 +0 17 

1 - 1U -0.14 

o.7518;1I 
0.4818:11 


1 9C+0.13 
1 - zo -0.12 

i.24i8;ll 

0-87l8;lo 
o.78i8;l 2 


2.89±8;g 
1.6818J? 

i-ool8:^ 
o.82i8;l 


< 15/0.19 (N) 

< 15/0.19 (E) 

> 15/0.19 (N) 

> 15/0.19 (E) 


U.Ol_ 24 
n 97+O.23 
u - z ' -0.20 

o.i8l8:8S 
o.o4l8:S 


1 01 +0.58 
i - ol -0.52 
1 47+O.59 
1 -*< -0.51 

0.3618.25 
0-82+0;!? 


o-67l8;i| 
0-80l8;i 9 5 

n 44+0.18 

u -^-0.14 

o-48l8:i 


29 +,u)y 

u - zy -0.08 
U - OU -0.09 

17+01° 

u -l ' -0.08 

14+ 016 

u - i *-0.13 


o.44l8;l^ 
o.58l8;l^ 
o.34l8;l 

n 49+0.40 

u '^ y -0.32 


49 +lUi 
u -^ y -o.n 

^7+o. n 

u -°'-0.09 

o.29i8;8? 

0.34±°;» 


1 26 +u - 68 

1 - zu -0.60 

1 3n +0 - 61 

1 - OU _0.47 

0.8818H 
1 20+ 154 

l- zu -0.96 



t Ratio of normalization Norm^ai/Norm^^. 



and 2T models give an Fe abundance to be ~ 0.7 solar 
and ~ 1.3 solar, respectively. These abundances decrease 
to about 0.3 solar at r ~ 0.3 ri 80 . In the outer region, the 
north and east fields give similar Fe abundances suggesting 
that the metals are mixed over a scale of 100 kpc , or 

~ 0.1 nso- 

The abundance profiles of the other elements are shown 
in Figure 4 (b)-(f). The O profile looks almost flat, al- 
though the statistical errors are fairly large. The IT model 
fit gives the O abundance to be ~ 0.3 solar, and the 2T 
model fit gives a two times that value. 

The abundance profiles of Mg, Si and S also show neg- 
ative gradients and resembles the Fe profile. All these 
abundances drop to about 1/2-1/3 of their central val- 
ues. The 2T model gives higher Mg and Si abundances in 
the central regions than the IT model, while the S abun- 
dance shows little dependence on the temperature model. 
Again, Mg, Si, and S abundances in the north and east 



fields are consistent with each other. 

The Ne and Ni abundances are significantly higher than 
those of other elements. Since the K-shcll lines of Ne and 
L-shell lines of Ni are completely mixed with the Fe-L 
lines, the derived abundances of Ne and Ni might include 
fairly large systematic uncertainties. 

Buote ct al. (2004) reported that the Fe abundance is 
about 0.15 solar at 0.2-0.4 r v ; r using the solar abundances 
by Grevesse & Sauval (1998), where the solar Fe abun- 
dance relative to H is 3.16xl0~ 5 by number. In contrast, 
we have derived about 0.3 solar at ~ 0.2 — 0.3 ri 80 based on 
the solar abundances of Lodders et al. (2003), where the 
solar Fe abundance relative to H is 2.95xl0 -5 by number. 
There is some difference in modeling the background (e.g., 
the Galactic emission), and the best-fit elemental abun- 
dance ratios in our analysis differs from that in Buote et 
al. (2004). This may cause the discrepancy in the Fe abun- 
dances. When the spectra at r > 15' were fitted with the 
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IT ICM model without the Galactic emission, we obtained 
the Fe abundance of 0.1-0.2 solar. However, this model 
were not able to reproduce the observed spectra at lower 
energy band, and \ 2 increased significantly. Another pos- 
sibility is that at the outer region of 0.3-0.4riso the Fe 
abundance might be very much lower than 0.3 solar. 

4-4- Profiles of abundance ratios of the ICM 

Figure 5 summarizes the abundance ratios of O, Mg, Si, 
and S divided by the Fe value, all taken in solar units. The 
abundance ratios, Mg/Fe, and S/Fe, are consistent with 
the absence of radial gradients. The O/Fe ratio is also 
consistent with having no radial gradient within rather 
large error bars. The discrepancies between the IT and 2T 
models become much smaller, except for the S/Fe ratio. 

The weighted average of the abundance ratios is calcu- 
lated for the central and offset regions and summarized 
in Table 4. On average, the abundance ratios are mostly 
consistent between the center and offset regions. The 2T 
model indicates O/Fe, Mg/Fe, Si/Fe, and S/Fe abundance 
ratios to be 0.6, 1.3, 0.7 and 1, respectively. The Ne/Fe 
ratio might have additional systematic uncertainty due to 
the overlap with Fe-L lines. 

4-5. Uncertainties in the spectral Fits 

To estimate the systematic errors, we varied the thick- 
ness of the OBF contaminant and the NXB levels by 
±10% in the same manner as in Sato et al. (2008b). The 
results are summarized in Figure 3 (a) and Figure 4. Here, 
the values in the north and east offset regions (r > 9') are 
averaged. Table 5 summarizes \ 2 values. The 10% de- 
crease in the OBF contaminant gives a slightly better \ 2 
value, while the change in the NXB level has little effect. 

We note that the systematic uncertainty in the derived 
temperature due to the change in the OBF contaminant is 
less than the statistical error in all the observed regions. In 
particular, the effect on the cool ICM component is almost 
negligible. The 10% change in the NXB level gives nearly 
the same ICM temperature and the effect is insignificant. 

The change in the metal abundance due to the 10% 
change in the OBF contaminant is smaller than 10-20% 
for all the elements, and it is very close to the statistical 
error. The change in the NXB level again has a negligible 
effect on the abundance determination. 

In summary, the uncertainty in the NXB level and the 
OBF contaminant have little effect on the temperature 
and abundance structures. 

4-6. Gas mass profile derived from the XMM-Newton ob- 
servation 

To derive an accurate gas mass profile, we used MOS 
and PN data from the XMM-Newton observation of 
NGC 5044. First, to derive accurate normalization for 
the hot and cool ICM components, we fitted the MOS and 
PN spectra, which are accumulated within thinner annuli, 
with the best-fit 2T model derived from the Suzaku spec- 
tral fit (3.1). We used the results of the central four an- 
nuli of the Suzaku spectra, or within <~ . 1 1 ri 80 ( 1 00 kpc) . 
The temperature and abundances of the ICM components 



were fixed to the best-fit values derived in Section 4.3, and 
normalizations of the vapec components were left free. To 
keep the spectral shape of the Galactic and the ICM emis- 
sion exactly the same as the Suzaku fit, we fixed the nor- 
malization ratio for the two components in the Galactic 
emission, as well as for the two ICM components, in all 
the annuli. 

We then calculated a volume filling factor f(r) of the 
cool component (Morita et al. 2006). The fractional vol- 
umes of the cool and hot gas, Vcooi and Vhot, to the total 
volume V are given by V coo i = fV and V hot = (1 - f)V, 
respectively. Because the normalizations of the 2T vapec 
models are expressed as Norm\ = C\i n\ V\ and Norrai = 
C12 n\ V2 using a certain common constant C12, the vol- 
ume filling factor f(r) can be solved under the pressure 
balance as (Morita et al. 2006), 



f(r) 



l+( Thot \ 2 N ° rm hot 

\T coo iJ Norm coo i 



(1) 



The radial profile of the filling factor is shown in Figure 
6 (left), for the XMM and Suzaku data. Considering 
the difference in the point spread function between the 
XMM and Suzaku data, the two results of the filling 
factor agree well. Within r < 0.07 rigo, the cool ICM 
component is more dominant than in the outer region. 
This filling factor profile has been fitted with a /3-modcl 
function, /(r) = [l + (r/r c j) 2 ] 3/3// ' 2 ; an d we obtained 
r cJ = 0.0055 ri 80 (4.85 kpc) and /?/ = 0.28. 

Based on the filling factor result, we derived the gas 
density distribution in the same way as in Morita et al. 
(2006) assuming that two temperature components of the 
ICM exist in the inner part of the extracted regions, or 
r <^ 0.11 ri 80 . In the outer region r >^ 0.11 ri 80 , it is 
not clear whether two ICM components exist. The derived 
radial profile of electron density n c is shown in Figure 6 
(right). The best-fit model for the filling factor is signifi- 
cantly smaller than the observed data within 0.05-0.2ri8o 
(Figure 6). However, this discrepancy does not affect the 
the total electron density within several %, since in this 
region, the electron density of the hot component domi- 
nates. Then, we fitted this density profile with a /3-modcl 

function, p gas (r) = Pgas^o [l + ( r / r c) 2 ] 3(3/2 , and obtained 
Pgas.o = 2.4 x 10 5 (M Q /kpc 3 ), r c = 0.032 r 180 (28 kpc) 
and = 0.53. 

The energy spectra for the XMM-Newton observation 
were also fitted with the same model as in our IT fitting 
(Sec. 3.1). As mentioned above, the single-temperature 
ICM component was assumed to exist in all extracted re- 
gions, or r <~ 0.17 riso- The temperature and abun- 
dances of the ICM component were fixed to the best-fit 
values derived in Section 4.3, and normalization of the 
vapec component was left free. We fitted the gas density 
profile, that was calculated from the normalizations de- 
rived from the IT spectral fit, with a /3-modcl function 
by fixing (3 = 0.51, which was obtained from the fitting of 
the XMM surface brightness profile. The derived param- 
eters of the /3-model are p ga s,o = 6-3 x 10 5 (M©/kpc 3 ) and 
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Fig. 5. Radial profiles of the abundance ratios of O, Mg, Si, and S divided by the Fe abundance (in units of solar ratio). The 
meanings of colors and symbols are the same as in Figure 4. 
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Fig. 6. (Left) Filling factor of the cool component derived from XMM (black) and Suzaku (gray), plotted along with the best-fit 
function of the results from XMM (black solid-line). The Suzaku data for the central four annuli, the north field, and east field are 
shown in solid, dashed and dotted lines, respectively. (Right) Electron density n c (cm -3 ) obtained using the results of the IT and 
2T fitting. The results for the cool (blue closed circle) and hot (red opened circle) ICM components and the total for the IT (gray) 
and 2T (black) models are plotted. The best-fit /3-model of the total 2T model is indicated by a solid black line. 
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Table 4. Weighted averages of the abundance ratios in solar units. 



Single temperature fitting 


field 


O/Fe 


Ne/Fe 


Mg/Fc 


Si/Fe 


S/Fe 


center 
offset 
all 


60 +lJ - iv 

u - ou -0.16 
77+0.36 
u - ' ' -0.35 

69+ 22 


1 fi9+ (J - 43 
L - oz -0A2 

2 Sfi+ 73 
i qq+o.42 


1 5fi+ (J - ib 

1 - oo -0.14, 

1 66+ 031 

l.DD_ 27 

1 61+ 017 


i.02lS:S 

0.78l°;£ 

o.9olS:8 


i - ou -0.15 
1 -"- )y -0.39 

1 60 +0 - 23 


Two temperature fitting 


field 


O/Fe 


Ne/Fe 


Mg/Fe 


Si/Fe 


S/Fe 


center 
offset 
all 


0.52^ 

o.6il8:l? 

0.561°;" 


i - oo -0.30 
o oo+0.85 
z.zz_ 7g 

9 no+0.46 
Z.UO_q 42 


1 - ZU -0.12 

1 42+ - 40 

1 -^ z -0.30 
i 01 +0.21 
1 - ol -0.16 


71 + U ' UY 

u - ' 1 -0.07 

61+° 18 
66+° 10 


S5+"- i,) 

u -°°-0.09 
-i i O+0.43 
1 - 1J -0.35 

0.99l° : 2 l 



Table 5. List of x 2 /dof for each fit to study systematic uncertainties 





Single 


temperature fitting 






Region 
(arcmin/ri 8 o) 


Nominal 


Contaminant 
+10% -10% 


NXB 
+ 10% -10% 


0-2/0-0.025 
2-4/0.025-0.050 
4-6/0.050-0.076 
6-9/0.076-0.11 


1574/899 
1441/899 
1194/899 
1027/899 


1685/899 
1514/899 
1239/899 
1055/899 


1483/899 
1389/899 
1160/899 
998/899 


1577/899 
1447/899 
1203/899 
1030/899 


1571/899 
1436/899 
1187/899 
1019/899 


< 15/0.19 (N) 

< 15/0.19 (E) 

> 15/0.19 (N) 

> 15/0.19 (E) 


1047/899 
963/899 
1091/899 
969/899 


1068/899 
999/899 
1095/899 
983/899 


1034/899 
937/899 
1094/899 
965/899 


1070/899 
990/899 
1100/899 
994/899 


1036/899 
949/899 
1095/899 
964/899 


Two temperature fitting 


Region 
(arcmin/ri 8 o) 


Nominal 


Contaminant 
+10% -10% 


NXB 
+10% -10% 


0-2/0-0.025 

2-4/0.025-0.050 

4-6/0.050-0.076 


1276/897 
1206/897 
1122/897 


1317/897 
1239/897 
1148/897 


1245/897 
1183/897 
1105/897 


1277/897 
1209/897 
1129/897 


1275/897 
1204/897 
1118/897 


6-9/0.076-0.11 


967/897 


981/897 


957/897 


973/897 


964/897 


< 15/0.19 (N) 

< 15/0.19 (E) 

> 15/0.19 (N) 

> 15/0.19 (E) 


999/897 
927/897 
1111/897 
967/897 


1001/897 
946/897 
1099/897 
975/897 


1002/897 
915/897 
1117/897 
965/897 


1008/897 
943/897 
1109/897 
986/897 


999/897 
921/897 
1115/897 
964/897 



r c = 0.017 mo (15 kpc). 

By integrating the gas density p g as, we derived the inte- 
grated gas mass profile for the IT and 2T models (Figure 
7) . For the 2T model, the sum of the cool and hot com- 
ponents is plotted. At r > 0.11 riso for the 2T fit and 
r > 0.17 rigo for the IT fit, we extrapolated the gas mass 
profile using the best-fit (3 model. A clear difference in 
the gas mass profile between the IT and 2T models is 
seen in the central region, r < ~ 0.03 ri 80 , while at 
r > ~ 0.03 riso, the two models give similar gas mass 
profiles. Within 2' from the center, the temperatures of 
the ICM derived from the IT model and from the cool 
component of the 2T model, which contributes most of 
the Fe-L emission, are nearly the same. In contrast, the 
IT model gives a lower Fe abundance than the 2T case. 
To reproduce similar intensity in the Fe-L band, a higher 
gas mass is required for the IT model. 

The gas mass profile derived from ROSAT observation 
(David ct al. 1994) is also shown in Figure 7. The ra- 



dial gradient of our gas mass profile agrees well with the 
ROSAT result, while the absolute values are larger than 
the ROSAT level by several tens of percentages, which 
is probably due to different assumptions about the metal 
abundance. 

Integrated and differential mass profiles of O, Mg, and 
Fe were derived from the gas mass and abundance profiles 
(Figure 7). In contrast to the gas mass profile, the metal 
mass profiles derived from the IT and 2T models agree 
well in the central region. Hereafter, we use the metal 
mass profiles derived from the 2T model. 

5. Discussion 

5.1. The Fe abundance profiles of the ICM 

The radial profile of Fe abundance in the NGC 5044 
group is compared with those in other groups, HCG 62 
(Tokoi et al. 2008), the NGC 507 group (Sato et al. 2008b), 
and the Fornax cluster (Matsushita et al. 2007a) observed 
with Suzaku (Figure 8). The three groups show similar 
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Fig. 7. (Left) Integrated mass of the hot gas (magenta), and those of O (red), Fe (blue), and Mg (green). Integrated B-band and 
K-band luminosities of galaxies, Lb and Lk (dotted lines), are also plotted. The dashed and solid lines correspond to the mass 
derived from the IT and 2T model fit, respectively. The gray dot-dashed line shows the gas mass derived using ROSAT (David et 
al. 1994). (Right) Differential metal mass density profiles (blue, red and green solid lines) derived for the 2T model and luminosity 
density profiles of Lb, and Lk- 



Fe abundance profiles. The Fe abundance is 1.2-1.5 so- 
lar at the center, and gradually decreases to 0.5 solar at 
O.lriso- The value of 0.5 solar is close to the level observed 
in general clusters of galaxies at 0.1-0.3ri 80 (Matsushita 
2008 ) . These groups are characterized by nearly symmet- 
ric X-ray morphology with a bright central galaxy, and 
the enhanced Fe abundance within 0.1 rigo should reflect 
Fe accumulation in the central galaxies over a long time. 
The NGC 5044 group shows similar Fe abundance in the 
north and east fields, suggesting that the group gas has 
evolved by keeping global spatial symmetry. In contrast, 
the Fornax cluster shows nearly constant Fe abundance 
at 0.5 solar over a radial range of 0.02-0.13 riso from the 
cD galaxy, NGC 1399. This cluster has asymmetric X-ray 
morphology, and NGC 1399 is significantly offset from the 
center (Paolillo et al. 2002, Scharf et al. 2005). Chandra 
observations suggest that NGC 1399 may be relatively 
moving against the ICM (Scharf et al. 2005). Therefore, 
around NGC 1399, a recent supply of metals from the cD 
galaxy produces a sharp peak of Fe abundance within 0.02 
riso- 

5.2. The abundance pattern of the ICM and contribu- 
tions from SN la and SN II 

Figure 9 summarizes the abundance pattern for O, Mg, 
Si, S and Fe in four systems: the NGC 5044 group, HCG 
62, the NGC 507 group and the Fornax cluster. The de- 
rived ratios, O/Fe, Mg/Fe, Si/Fe and S/Fe, all scatter 
around unity, i.e., the solar ratio. The four systems show 
similar Si/Fe and S/Fe ratios. The NGC 5044 group shows 
the smallest O/Fe ratio, 0.6, compared with ~ 0.8 in the 
other three systems. In contrast, the NGC 5044 group 
shows the highest Mg/Fe ratio of 1.2 among the four sys- 
tems derived from the 2T model. The 2T and IT models 
give different abundance ratios by a few tens of percent- 
ages. The uncertainty due to coupling with the Fe-L lines 
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Fig. 8. Radial profiles of the Fe abundance in the NGC 
5044 group (diamonds), the Fornax cluster (Matsushita et al. 
2007a, open circles), the NGC 507 group (Sato et al. 2008b, 
closed circles), and HCG 62(Tokoi et al. 2008, closed trian- 
gles). 

may result in a 20-30% systematic error in the Mg abun- 
dance (Matsushita et al. 2007a), and the estimation of the 
Galactic emission leaves some uncertainty in the O abun- 
dance. Considering the systematic uncertainty mentioned 
here, the abundance pattern for O, Mg, Si, S, and Fe does 
not differ significantly among the four poor systems. 

As shown by Sato et al. (2007b), the abundance pattern 
of the NGC 5044 ICM within 0.1ri 80 and 0.1-0.3ri 80 can 
be fit by a combination of average SNe la and SNc II 
yields (Figure 10). The SNe la and II yields were taken 
from Iwamoto et al. (1999) and Nomoto et al. (2006), 
respectively. About 80% of Fe and ~ 40% of Si and S in 
the ICM are synthesized by SN la. The derived number 
ratio of SNe II to SNe la is 3.1 ± 0.6 and 3.3±0.7, for r < 
O.lrigo and O.lrigo < r < 0.3rigo, respectively, assuming 
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Fig. 9. Weighted averages of the abundance ratios of vari- 
ous elemental species compared to Fe for the central (red; 
r <~ O.lriso) and offset (blue; r >~ O.lriso) fields of the NGC 
5044 group. Those in the regions within ~ O.lriso (magenta) 
and 0.1 ~ 0.3ri8o (light blue) of the NGC 507 group (Sato et 
al. 2008b), HCG 62 (Tokoi et al. 2008), the Fornax cluster 
(Matsushita et al. 2007a), and an elliptical galaxy NGC 1404 
(Matsushita et al. 2007a) are also plotted. 

a classical deflagration model, W7, for SNe la. This is 
close to the value of 3-4, obtained for the other groups of 
galaxies, NGC 507 and HCG 62, and several other clusters 
of galaxies (Sato et al. 2007b). 

When systems have a similar initial mass function of 
stars (IMF) and are old enough, most SN la and SN II 
would have already exploded. To explain abundance pat- 
tern of stars in the solar neighborhood, lifetimes of SN 
la are confined within 0.5-3 Gyr, with typical lifetime of 
1.5 Gyr (Yoshii et al. 1996). In clusters of galaxies, to 
account for Fe mass in the ICM, the past average rate 
of SNe la was much larger than the present rate in el- 
liptical galaxies (e.g., Renzini et al. 1993). This result 
indicates that lifetimes of most of SN la are much shorter 
than the Hubble time. In this case, the final abundance 
pattern should be similar. If the IMF is close to that in 
our Galaxy and most of stars in our Galaxy and clusters 
were already formed before a few Gyrs ago, the abundance 
pattern should naturally be similar to the solar abundance 
pattern. Within O.lriso of NGC 5044, the gas mass and 
stellar mass are pretty close (§5.3). The mean stellar 
metallicity and [a/Fe] derived from optical observations 
are about 1 solar and 0.34 ±0.17, respectively (Kobayashi 
& Arimoto 1999, Annibali et al. 2007). Therefore, the 
sum of the low O /Fe ratio found in the ICM and the high 
stellar a/Fe values may give roughly the solar abundance 
pattern. However, in this study, we are combining the 
abundances of two distinct bodies, stars and ISM, which 
may have very different enrichment histories. 

5.3. The radial profiles of the metal mass-to-light ratios 

Since metals in the ICM are all synthesized in galaxies, 
the metal-mass-to-light ratio is a useful measure for the 
study of ICM chemical evolution. To estimate the metal 
mass-to-light ratio, we calculated B-band and K-band lu- 



minosity profiles. Ferguson & Sandage (1990) (hereafter, 
FS90) have optically studied the NGC 5044 group, and 
identified 162 member galaxies whose apparent magni- 
tudes, me, are brighter than me ~ 20, with a complete- 
ness limit of tob ~ 18. They also show that the central 
elliptical galaxy NGC 5044 has an apparent magnitude 
of me = 11.9, or logLe/iB,© = 10.7 using the luminos- 
ity distance Dl = 38.9 Mpc and the foreground Galactic 
extinction A B = 0.300 (Schlegel et al. 1998) from NED. 
Since this galaxy dominates the central region, we ap- 
plied the de Vaucouleurs' law to calculate the luminosity 
profile within 0.076 riso (6'), and deprojected the profile 
using the expression by Mellier & Mathez (1987). In the 
outer region, r > 0.076 r^o, we integrated the luminosi- 
ties of the member galaxies from FS90 and deprojected 
the profile. 

Since the K-band luminosity of a galaxy correlates well 
with the stellar mass, we also collected K-band mag- 
nitudes of galaxies in a box of 2 x 2 deg 2 centered on 
NGC 5044 from the Two Micron All Sky Survey (2MASS). 
NGC 5044 itself has an apparent magnitude bik = 7.845, 
or logl/K/^K,© = 11-4 using the foreground Galactic ex- 
tinction A K = 0.026 (Schlegel et al. 1998) from NED. The 
average surface brightness in the region of r > riso = 79.2', 
i'K.bgd = 3.01 x 10 9 Lk.q, is subtracted as the background. 
We then calculated the K-band luminosity profile in the 
same way as for the B-band profile. 

Forty-three objects are detected in both the B-band and 
K-band (Figure 11). As shown in Figure 11, the two mag- 
nitudes of these objects correlate well with each other. 

The integrated B-band and K-band luminosity profiles 
are shown in the left panel of Figure 7 and in Figure 12. 
The two band profiles resemble each other; the ratio of 
the integrated luminosity, Lb/Lk, is nearly constant at 
^0.22 (figure 12). This corresponds to B — K ^ 4.1, which 
is consistent with B — K = 4.2 for early- type galaxies in 
Lin & Mohr (2004). The NGC 5044 galaxy dominates the 
luminosity profile in the region r <^ 0.07 riso, while in 
the outer region, r >^ 0.07 riso, the contribution from 
fainter galaxies in the group becomes important. 

The luminosity density profiles of Lb and Lk are plot- 
ted in the right panel of Figure 7. These profiles are com- 
pared with the mass density profiles of O, Mg and Fe from 
the 2T model fit. We note that the mass densities of stars 
and O both indicate a shoulder-like structure in the region 
from 0.06riso to 0.2ri 8 o- The profiles of Mg and Fe do not 
show such a shoulder and decrease smoothly. Therefore, 
the oxygen distribution seems to closely correlate with 
that of small galaxies in the group. 

The integrated mass-to-light ratios for O, Mg and Fe 
(OMLR, MMLR and IMLR) using B-band and K-band lu- 
minosities are summarized in Table 6 and Figure 13 (left). 
The error bars of the mass-to-light ratios include only the 
abundance errors. The profiles increase with the radius up 
to r ~ 0.1 riso, in the region where NGC 5044 dominates, 
and then become almost constant in the outer region. 

Our IMLR value with the B-band luminosity is ~ 
3.6 x 10~ 3 expressed in terms of solar values at r < 0.3 riso- 
This value is consistent with the result of Buote et al. 
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Fig. 10. Fit results of number ratios of elements to Fe of the NGC 5044 group within 0.1 rigo (left) and 0.1—0.3 rigo (right). 
Top panel shows the number ratios (closed circles). Blue dashed and dotted lines correspond to the contributions of SNe la (W7) 
and SNe II, respectively, and the red line corresponds to the sum of the two. Ne (atomic number = 10) is excluded from the fit. 
The bottom panel indicates fractions of the SNe la contribution to total mass in the best-fit model for each element. 
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Fig. 11. (Left) Objects detected in B-band and K-band (black closed circle: B-band, red open circle: K-band). (Right) 
Correlation of magnitude of B-band and K-band for objects detected in both bands. 
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Fig. 12. Integrated luminosity profiles in B-band and 
K-band. The dashed lines indicate luminosity of the central 
galaxy NGC 5044. 

(2004), since our larger derived Fe abundance gives a 
smaller gas mass. 

We also calculated differential metal-mass-to-light ra- 
tios, as shown in the right panel of Figure 13. The dif- 
ferential OMLR becomes almost flat in the region r >^ 
0.07 ri8o, while IMLR and MMLR show a clear drop in 
0.1-0.2riso compared with the level at 0.07ri 80 . In other 
words, the O distribution seems to be more extended than 
those of Fe and Mg, although the error in the data is fairly 
large. There may be some difference in the metal pro- 
duction process between the giant elliptical galaxy, NGC 
5044, and smaller galaxies, in connection with a different 
history of SN la and SN II activity. 

5.4- Comparison of the metal-mass-to-light ratios with 
other systems 

The metal-mass-to-light ratios of the NGC 5044 group 
are compared with those in other groups and clusters. 
The objects are, the NGC 507 group and HCG 62, the 
Fornax cluster as a poor cluster, and other clusters of 
galaxies with ICM temperatures of 2-4 keV: A 262, A 
1060, AWM7, and the Centaurus cluster, as summarized 
in Figure 14 and Table 7. The gas mass of these systems 
were derived using XMM-Newton data except the Fornax 
cluster, where the ROSAT data was used. 

The OMLR, MMLR, and IMLR for clusters of galaxies 
with kT ~ 2 — 4 keV tend to be higher than those in groups 
of galaxies and the poor cluster of galaxies with kT ~ 1 
keV. However, the data points might have some system- 
atic uncertainty, possibly because of different observations 
of £?-band luminosity of galaxies, statistics in the number 
counts of galaxies, and systematic shift of gas mass de- 
rived from XMM observations with different sensitivities 
and fields of view. 

The difference in the OMLR values between groups and 
clusters is about a factor of 3-6, and systematically larger 
than the IMLR and MMLR difference, which is a factor 
of 2-3. Among groups, the NGC 5044 group and HCG 62 



have the highest IMLRs, which are comparable to those 
of clusters. In contrast, OMLR of the NGC 5044 appears 
to be lower than those of clusters. Among the groups and 
poor clusters, the Fornax cluster shows the lowest IMLR 
and OMLR. In the Fornax cluster, the asymmetric X- 
ray emission may point to large-scale dynamical evolution, 
which might have hampered the strong concentration of 
hot gas in the center. 

These poor systems also differ from richer systems in 
that the gas density profile in the central region is signifi- 
cantly flatter and the relative entropy level is correspond- 
ingly higher (e.g. Ponman et al. 1999). These peculiar 
central characteristics in the poor systems are thought to 
be best attributed to an injection of energy (preheating) 
into the gas before clusters have collapsed (e.g. Kaiser 
1991). Note also that more recent observations and sim- 
ulations of entropy profiles indicate some modification to 
this simple picture (e.g. Ponman et al. 2003). 

Metal distribution in the ICM can be a powerful tracer 
of the history of such a gas heating in the early epoch, 
since the relative timing of metal enrichment and heat- 
ing should affect the present amount and distribution of 
metals in the ICM. To look into this, we assume that all 
galaxies synthesize a similar amount of metals per unit 
stellar mass. If metal enrichment occurred before the en- 
ergy injection, the poor systems would carry relatively 
smaller metal mass with smaller gas mass than those in 
rich clusters, while metal abundance would be quite sim- 
ilar to those in rich clusters. In contrast, if metal enrich- 
ment occurred after the energy injection, the metal mass 
becomes comparable to those in rich clusters and indicates 
a higher abundance. The difference in O and Fe distribu- 
tions would constrain the timing of energy injection, since 
O was synthesized by SN II during galaxy star formation, 
while Fe was mainly produced by SN la much later than 
this period. 

The OMLR shows a larger difference in clusters and 
groups than the IMLR does, which may reflect the con- 
dition that SN II products are more extended in groups 
due to the effect of preheating. Fe was synthesized later 
by SN la, and may be irrelevant to preheating. To fur- 
ther investigate the problem of heating and enrichment, 
we need numerical simulations including the effect of the 
preheating and metal enrichment by both SN II and SN 
la. 

6. Summary and Conclusion 

Suzaku observations of the NGC 5044 group, with a 
relatively symmetric X-ray morphology, allowed determi- 
nation of the abundances of O, Mg, Si, S and Fe in the 
ICM up to 0.3ri8o fairly accurately. The Fe abundance 
around NGC 5044 is about a solar and drops to 0.3 solar 
at 0.3ri 80 . The abundance ratios, Mg/Fe, Si/Fe and S/Fe, 
are close to the solar ratio, while O/Fe is 0.5-0.6 in solar 
units. The abundance pattern indicates that most of the 
Fe should have been synthesized by SN la. The O distri- 
bution closely resembles the light distribution of galaxies, 
in particular in the outer region. The OMLR of the NGC 
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Table 6. Summary of metal-mass-to-light ratios of NGC 5044 
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Fig. 13. (Left) Radial profiles of integrated IMLR (blue), OMLR (red) and MMLR (green) in B-band (solid lines) and K-band 
(dashed lines). (Right) Radial profiles of differential mass-to-light ratios with the same representation as in the left panel. 



5044 group is lower than those of clusters, but the IMLR 
is comparable to those of clusters. The metal distribution 
of the ICM can be used as a tracer of the past history of 
heating and enrichment in the clusters of galaxies. 
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